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Abstract 

We study the qualitative properties of a limiting elliptic system arising in phase separation 
for Bose-Einstein condensates with multiple states: 




When n = 1, we prove uniqueness of the one-dimensional profile. In dimension 2, we prove 
that stable solutions with linear growth must be one-dimensional. Then we construct entire 
solutions in K 2 with polynomial growth \x\ d for any positive integer d > 1. For d > 2, these 
solutions are not one-dimensional. The construction is also extended to multi-component 
elliptic systems. 
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1 Introduction and Main Results 

Consider the following two-component Gross-Pitaevskii system 

- Au + an 3 + Av 2 u = Xiu in O, (1.1) 
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- Av + f3v 3 + Au 2 v = X 2 v in 0, (1.2) 

u > 0, u > in O, (1.3) 

u = 0, v = ondn, (1.4) 



^ = iVi, / ^ = 7V 2 , (1.5) 

where a,/3,A > and £1 is a bounded smooth domain in M. n . Solutions of (1. 1)-(1.5) can be 
regarded as critical points of the energy functional 



Ea(u,v) 



[ (\Vu\ 2 + \Vv\ 2 ) + V + ^4 + A 2 u 2 j (16) 
Jn 2 2 2 

on the space (u,v) £ Hq(Q) x Hq(Q) with constraints 

u 2 dx = N l , t v 2 dx = N 2 . (1.7) 



The eigenvalues Aj's are Lagrange multipliers with respect to (1.7). Both eigenvalues Xj = 
Aj,A>i = 1) 2, and eigenfunctions u = u\, v = v\ depend on the parameter A. As the parameter A 
tends to infinity, the two components tend to separate their supports. In order to investigate the 
basic rules of phase separations in this system one needs to understand the asymptotic behavior 
of (u\,v\) as A — > +oo. 

We shall assume that the solutions (u\,v\) of (1.1)-(1.5) are such that the associated eigen- 
values A^a's are uniformly bounded, together with their energies E\(u\,v\). Then, as A — >■ +oo, 
there is weak convergence (up to a subsequence) to a limiting profile (uoo^oo) which formally 
satisfies 

J -Aiioo + m4, = A^ooUoo in tt u , ^ 
| - Avoo + fivlc = A 2 ,oo'Voo in Q v , 

where Q u = {x E Q : Uoo(x) > 0} and 1!„ = {i £ fl : Voo(x) > 0} are positivity domains composed 
of finitely disjoint components with positive Lebesgue measure, and each Aj j00 is the limit of A^a's 
as A — > oo (up to a subsequence) . 

There is a large literature about this type of questions. Effective numerical simulations for (1.8) 
can be found in [5], [6] and [13]. Chang-Lin-Lin-Lin [13] proved pointwise convergence of (u\,v\) 
away from the interface r = {x € : u OD (x) = Voo(x) = 0}. In Wei-Weth [27] the uniform 
equicontinuity of (ua, v\) is established, while Noris-Tavares-Terracini-Verzini [24] proved the 
uniform-in- A Holder continuity of (u\,v\). The regularity of the nodal set of the limiting profile 
has been investigated in [12, 26] and in [16]: it turns out that the limiting pair (u 00 (x),v 00 (x)) is 
the positive and negative pair (w + , w~) of a solution of the equation —Aw + a(w + ) 3 — [3(w~)' i = 

Al,ooW + - \2,odW~. 

To derive the asymptotic behavior of (u\,v\) near the interface T = {x £ £1 : Uoo(x) = 
Voo(x) = 0}, one is led to considering the points x\ G f2 such that ua(xa) = v\(x\) = m\ — > 
and xa — > € 7 C £1 as A — > +00 (up to a subsequence). Assuming that 

m\k -> C > 0, (1.9) 
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(without loss of generality we may assume that Cq = 1), then, by blowing up, we find the following 
nonlinear elliptic system 

Au = uv 2 , Av = vu 2 , u,v>0 in R n . (1.10) 

Problem (1.10) has been studied in Berestycki-Lin- Wei-Zhao [8], and Noris-Tavares-Terracini- 
Verzini [24]. It has been proved in [8] that, in the one-dimensional case, (1.9) always holds. 
In addition, the authors showed the existence, symmetry and nondegeneracy of the solution to 
one-dimensional limiting system 

" 2 " 2 

u =uv ,v = vu , u, v > in M. (1-11) 

In particular they showed that entire solutions are reflectionally symmetric, i.e., there exists 
xo such that u(x — Xo) = v(xq — x). They also established a two-dimensional version of the De 
Giorgi Conjecture in this framework. Namely, under the growth condition 

u(x) + v(x) <C(l + \x\), (1.12) 

all monotone solution is one dimensional. 

On the other hand, in [24], it was proved that the linear growth is the lowest possible for 
solutions to (1.10). In other words, if there exists a € (0, 1) such that 

u(x) +v(x) < C(l + \x\) a , (1.13) 

then u, v = 0. 

In this paper we address three problems left open in [8]. First, we prove the uniqueness of 
(1.11) (up to translations and scaling). This answers the question stated in Remark 1.4 of [8]. 
Second, we prove that the De Giorgi conjecture still holds in the two dimensional case, when we 
replace the monotonicity assumption by the stability condition. A third open question of (1.10) 
is whether all solutions to (1.10) necessarily satisfy the growth bound (1.12). We shall answer 
this question negatively in this paper. 

We first study the one-dimensional problem (1.11). Observe that problem (1.11) is invariant 
under the translations (u(x),v(x)) — > (u(x + t),v(x + t)), Vi 6 R and scalings (u(x),v(x)) — > 
(Att(Ax), Av(\x)), VA > 0. The following theorem classifies all entire solutions to (1.11). 

Theorem 1.1. The solution to (1.11) is unique, up to translations and scaling. 

Next, we want to classify the stable solutions in M 2 . We recall that a stable solution (u,v) to 
(1.10) is such that the linearization is weakly positive definite. That is, it satisfies 

/ [V<£| 2 + |VVf + v 2 ip 2 + uV + 4uv<pi>] > 0, V^i e c °°(R n ). 

In [8], it was proved that the one-dimensional solution is stable in M n . Our first result states 
that the only stable solution in M 2 , among those growing at most linearly, is the one-dimensional 
family. 
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Theorem 1.2. Let (u,v) be a stable solution to (1-10) in R 2 . Furthermore, we assume that the 
growth bound (1-12) holds. Then (u,v) is one- dimensional, i.e., there exists a € M 2 , \a\ = 1 such 
that (u, v ) = (U(a ■ x), V(a ■ x)) where (U, V) are functions of one variable and satisfies (1.11 ). 

Our third result shows that there are solutions to (1-10) with polynomial growth \x\ d that 
are not one dimensional. The construction depends on the following harmonic polynomial <3? of 
degree d: 

$ := Re(z d ). 

Note that <3? has some dihedral symmetry; indeed, let us take its d nodal lines Li, ■ ■■ ,Ld and 
denote the corresponding reflection with respect to these lines by T%, • • • ,7^. Then there holds 

$(TiZ) = (1.14) 

The third result of this paper is the following one. 

Theorem 1.3. For each positive integer d > 1, there exists a solution (u,v) to problem (1.10), 
satisfying 

1. u-v > in {$ > 0} and u - v < in {$ < 0}; 

2. u>& + and v > 

3. Vi = 1, ■ ■ ■ ,d, u(Tiz) = v(z); 

4- Vr > 0, the Almgren frequency function satisfies 

r fn m\ |Vu| 2 + \Vv\ 2 + U 2 V 2 

N(r) := JML. -^— 2 < d; (1.15) 

JdB r (0) U V 

5. 

lim N(r) = d. (1.16) 

Note that the one-dimensional solution constructed in [8] can be viewed as corresponding to 
the case d = 1. For d > 2, the solutions of Theorem 1.3 will be obtained by a minimization 
argument under symmetric variations (<p,ip) (i.e. satisfying (p o Tj = ip for every reflection Tj). 
The first four claims will be derived from the construction. See Theorem 4.1. 

Regarding the claim 5, we note that by Almgren's monotonicity formula, (see Proposition 
5.2 below), the Almgren frequency quotient N(r) is increasing in r. Hence lim r ^ +00 N(r) exists. 
To understand the asymptotics at infinity of the solutions, one way is to study the blow-down 
sequence defined by: 

1 1 

{u R (x),v R (x)) := (T^y u ( jRx )7^y ?; (- Rx ))' 

where L(R) is chosen so that 



Idi 



'SBi(O) 

In Section 6, we will prove 



Ur + Vr = 1. 
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Theorem 1.4. Let (u,v) be a solution of (1.10) such that 



d := lim N(r) < +00. 

Then d is a positive integer. As R — >• 00, (ur,vr) defined above (up to a subsequence) converges 
to (^ + , uniformly on any compact set of~R N where ^ is a homogeneous harmonic polynomial 
of degree d. If d = 1 then (it, v) is asymptotically flat at infinity. 

In particular this applies to the solutions found by Theorem 1.3 to yield the following property 

Corollary 1.5. Let (u,v) be a solution of (1.10) given by Theorem 1.3. Then 

(u R (x),v R (x)) := (^u(Rx)^v(Rx)) 

converges uniformly on compact subsets ofM 2 to a multiple 0/ (<1? + , <I> _ ), where <5 := Re(z d ). 

Theorem 1.4 roughly says that (u,v) is asymptotic to (\E ,+ , x I /_ ) at infinity for some homo- 
geneous harmonic polynomial. The extra information we have in the setting of Theorem 1.3 is 
that \& = <J? = Ke(z d ). This can be inferred from the symmetries of the solution (property 3 in 
Theorem 1.3). 

For another elliptic system with a similar form, 

f Ait = uv.u > in R n , 

< (1.17) 
[Au = TO,o0inR n 

the same result has been proved by Conti-Terracini-Verzini in [15]. In fact, their result hold for any 
dimension n > 1 and any harmonic polynomial function on M n . Note however that the problem 
here is different from (1.17). Actually, equation (1.17) can be reduced to a single equation: indeed, 
the difference u — v is a harmonic function (A (it — v) = 0) and thus we can write v = u — $ where 
$ is a harmonic function. By restricting to certain symmetry classes, then (1.17) can be solved 
by sub-super solution method. However, this reduction does not work for system (1.10) that we 
study here. 

For the proof of Theorem 1.3, we first construct solutions to (1.10) in any bounded ball -Br(O) 
satisfying appropriate boundary conditions: 

An = ni; 2 , in B R (0), 

Av = vu 2 , in Bfl(0), (1.18) 
u = <$ + ,v = <f>~ ondB R (0). 

This is done by variational method and using heat flow. The next natural step is to let 
R — > +00 and obtain some convergence result. This requires some uniform (in R) upper bound 
for solutions to (1.18). In order to prove this, we will exploit a new monotonicity formula for 
symmetric functions (Proposition 5.7). We also need to exclude the possibility of degeneracy, 
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that is that the limit could be or a solution with lower degree such as a one dimensional 
solution. To this end, we will give some lower bound using the Almgren monotonicity formula. 

Lastly, we observe that the same construction works also for a system with many components. 
Let d be an integer or a half-integer and 2d = hk be a multiple of the number of components k, 
and G denote the rotation of order 2d. In this way we prove the following result 

Theorem 1.6. There exists a positive solution to the system 

k 

Aui = Ui uj, in C = R 2 , i = 1, . . . , k, 

Ui > 0, i = 1, . . . , k, 
having the following symmetries (here ~z is the complex conjugate of z) 



(1.19) 



Furthermore, 



and 



Ui(z) = Ui(G h z), on C , i = 1, . . . , k, 

Ui{z) = m+i(Gz), on C , % = 1, . . . , k, 

u k+ i(z) = ui(z), on C 

Uk+2-i(z) = Ui(z), onC,i = l,...,k. 

If k 

l ™ 1T2d / Yl u i = b G (°> +°°) ' 



(1.20) 



r/ Br(0) E?|v^| 2 + Ei<^H A 

hm — : = d . 



JdB r (0) Z^l U i 



ldB r (0) 

The problem of the full classification of solutions to (1.10) is largely open. In view of our 
results, one can formulate several open questions. 

Open problem 1. We recall from [8] that it is still an open problem to know in which dimension 
it is true that all monotone solution is one-dimensional. A similar open question is in which 
dimension it is true that all stable solution is one-dimensional. We refer to [2], [20], [18], [23], and 
[25] for results of this kind for Allen-Calm equation. 

Open problem 2. Let us recall that in one space dimension, there exists a unique solution to 
(1.11) (up to translations and scalings). Such solutions have linear growth at infinity and, in the 
Almgren monotonicity formula, they satisfy 

lim N(r) = 1. (1.21) 

r— >-+oo 

It is natural to conjecture that, in any space dimension, a solution of (1.10) satisfying (1.21) is 
actually one dimensional, that is, there is a unit vector a such that (u(x), v(x)) = (U(a-x), V(a-x)) 
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for x £ W l , where (U, V) solves (1.11). However this result seems to be difficult to obtain at this 
stage. 

Open problem 3. A further step would be to prove uniqueness of the (family of) solutions 
having polynomial asymptotics given by Theorem 1.3 in two space dimension. A more challenging 
question is to classify all solutions with 

lim N(r) = d. (1.22) 



Open problem 4. For the Allen-Cahn equation Au + u — u = in R 2 , solutions similar to 
Theorem 1.3 was first constructed in [17] for d = 2 and in [1] for d > 3. (However all solutions to 
Allen-Cahn equation are bounded.) On the other hand, it was also proved in [19] that Allen-Can 
equation in R 2 admits solutions with multiple fronts. An open question is whether similar result 
holds for (1.10). Namely, are there solutions to (1.10) such that the set {u = v} contains disjoint 
multiple curves? 

Open problem 5. This question is related to extension of Theorem 1.3 to higher dimensions. 
We recall that for the Allen-Cahn equation Am + u — u = in M 2m with m > 2, saddle-like 
solutions were constructed in [10] by employing properties of Simons cone. Stable solutions to 
Allen-Cahn equation in M 8 with non planar level set were found in [23], using minimal cones. We 
conjecture that all these results should have analogues for (1.10). 



2 Uniqueness of solutions in R: Proof of Theorem 1.1 

In this section we prove Theorem 1.1. Without loss of generality, we assume that 

lim u{x) = +oo, lim v(x) = 0. (2-1) 

The existence of such entire solutions has been proved in [8]. By symmetry property of solutions 
to (1.11) (Theorem 1.3 of [8]), we may consider the following problem 

{II n II n 

u = uv , v = vu , u, v > in R, 
(2.2) 
lim u (x) = — lim v (x) = a 
x— >+oo x-^—oo 

where a > is a constant. We now prove that there exists a unique solution (u, v) to (2.2), up to 
translations. We will prove it using the method of moving planes. 

First we observe that for any solution (u, v) of (2.2), u and v decay exponentially at infinity. 
Integration shows that as x — > +oo, \u (x) — a\ decays exponentially. (See also [8].) This implies 
the existence of a positive constant A such that 

\u(x) — ax + \ + \v(x) — ax~\ < A. (2-3) 
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Moreover, the limits 

lim (u(x) — ax + ), lim (v(x) — ax~) 

X— ¥+00 X— > — oo 

exist. 

Now assume (u\,vx) and (u 2 ,v 2 ) are two solutions of (2.2). For t > 0, denote 

:= + t),vi )t (x) := v\(x + t). 

We want to prove that there exists an optimal to such that for all t > to, 

u\,t{x) > U2(x),vi 7 t(x) < v 2 (x) in R. (2.4) 

Then we will show that when t = to these inequalities are identities. This will imply the uniqueness 
result. 

Without loss of generality, assume (in, i>i) and (7/2,^2) satisfy the estimate (2.3) with the same 
constant A. 

Step 1. For t > ^ (A as in (2.3)), (2.4) holds. 

Firstly, in the region {x > —t + ^r}, by (2.3) we have 

ui,t(x) >a(x + t) - A> ax + + A > u 2 (x); (2.5) 
while in the region {x < — t + — }, we have 

vi,t(x) < a(x + t)~ + A < ax~ - A < v 2 (x). (2.6) 
On the interval {x < — t + — }, we have 



u i,t = u i,tvf )t < u i,tvl, 

" 2 
u 2 = u 2 v 2 . 



With the right boundary conditions 



2A 2A 

ui,t(—t-\ )>u 2 (-t-\ ), lim ui t (x)= lim u 2 (x) = 0, 

a a x—>—co x— >— 00 

a direct application of the maximum principle implies 

inf (ux,t - u 2 ) > 0. 
{ ,<-i+M } 

By the same type of argument also show that 



sup (vi >t - v 2 ) < 0. 
{x>-t+m 



(2.7) 



Therefore, we have shown that for t > ^^,«it > u 2 and v\j < v 2 . 



Step 2. We now decrease the t to an optimal value when (2.4) holds 

t = inf {t | such that (2.4) holds for all t > t'}. 

Thus to is well defined by Step 1. Since — (tii,t„ — Ut i)' + v xt ( u i,to ~ u z) — Oj ~( v 2 — v i,t )" + 
u 1t ( v 2 ~ v i,to) — 0) by the strong maximum principle, either 

Ui,to(x) = u 2 (x),v 1)to (x) = v 2 (x) in R, 



or 



ui,t (x) > u 2 (x),v ltto (x) < v 2 (x) in R. (2.J 



Let us argue by contradiction that (2.8) holds. By the definition of to, there exists a sequence 

sither 

inf - u 2 ) < 0, (2.9) 



of t k < to such that lim t k = to and either 



or 

sup(«i t , - v 2 ) > 0. 

K 

Let us only consider the first case. 

Define u>i tk '■= Uij k — u 2 and w 2>k '■= v 2 — i>i,t fe - Direct calculations show that they satisfy 



w i,k + v i,t k w i,k = u 2 (v 2 + v lttk )w 2tk in 



I - w 2,k + u i,t h w 2,k = v 2 (u 2 + u ljtk )w 1)k in R. 
We use the auxiliary function g(x) = log(|x| + 3) as in [14]. Note that 

g > 1, g" <0 in {x ± 0}. 
Define Wi jk := w\^jg and w 2ik := w 2jk /g. For i/Owe have 



(2.10) 



(2.11) 



-w 2k - 2—w 2 k + [u\ t ]w 2k = v 2 (u 2 + Uij k )wi,k, in R- 

By definition, it?i^ and u^.fc are bounded in R, and hence 

^l,fe) w 2 ^k — > as | x | — ^ oo . 

In particular, in view of (2.9), we know that 'vcd-^{wi tk ) < is attained at some point x k ,\. 

Note that \x k ,i\ must be unbounded, for if x k> i — > Xoo,t k — > to, then tui^a^i) — > iti,t (xoo) 
U2(xoo) = 0. But this violates the assumption (2.8). 
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Since \xk,i\ is unbounded, at x = Xk,i there holds 

w 1 k > and w lk = 0. 
Substituting this into the first equation of (2.11), we get 

[vi,t k (xk.i) 2 ~ 9 , Xh '\ jwi^ixk ,i) > u 2 (x kA )(v 2 (xk,i) + vi t h {x k l ))w 2! k(xk,i) (2.12) 

9{ x k,l) 

which implies that w 2 k(xk i) < 0. Thus we also have inf w 2 k < 0. Assume it is attained at xu 2 . 

' R 

Same argument as before shows that \x k ,2\ must also be unbounded. Similar to (2.12), we have 
g"r Xk 2 \ 

[u\,t k {x k , 2 ) 2 -, — V]^2,fc(2;fc,2) > v 2 {x k2 ) (u 2 (x k , 2 ) + u lttk (x k)2 ))wik(xk, 2 ). (2.13) 

Observe that 

w 2;k (x k , 2 ) = inf w 2 , k < w 2 , k (x k>1 ), 

K 

wik{x k l ) = inf w l k < w lk (x k2 ). 
Substituting these into (2.12) and (2.13), we obtain 

~ / \ . U 2 {x k ,l)[v2{Xk,l) + Vl,t k (x k ,l)] V2(Xk,2)[u 2 (x k , 2 ) + U htk (x kt2 )] 

wi,k{ x k,i) > 7T, ; ; Wi,ki x k,V- ( 2 - 14 ) 

\2 9 (Sfc.i) „ („ \2 9 (Xk,2) 



Since Wi >k (x kj i) < 0, we conclude from (2.14) that 

> 1 (2.15) 



"2(^,1)^2(^,1) + V l >tk (X k> i)] V 2 (x k>2 )[u 2 {x kt2 ) + U 1)tk (x k ,2)] 



II. 

where \x k ±\ -> +00, |x fc)2 | -> +00. This is impossible since ~ - upwurp^ as M -> +00, 

and we also use the decaying as well as the linear growth properties of u and t> at 00. 

We have thus reached a contradiction, and the proof of Theorem 1.1 is thereby completed. 



3 Stable solutions: Proof of Theorem 1.2 

In this section, we prove Theorem 1.2. The proof follows an idea from Berestycki-Caffarelli- 
Nirenberg [7]-see also Ambrosio-Cabre [2] and Ghoussoub-Gui [20]. First, by the stability, we 
have the following 

Lemma 3.1. There exist a constant A > and two functions (p > and ip < 0, smoothly defined 
in M 2 such that 

A(p = v 2 (p + 2uvtf) — \(p, 

2 (^*1) 

= 2uv(p + v if) — Xijj. 
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Proof. For any R < +00 the stability assumption reads 

Jr rm l V( ^l 2 + I W| 2 + v 2 p 2 + w 2 ^ 2 + 4uv(f»f} 
X(R) := min — ^ = = > 0. 

V>,4'&H^B R (0))\{0} J Br{0) p 2 + V 

It's well known that the corresponding minimizer is the first eigenfunction. That is, let ((pR,ip R ) 
realizing X(R), then 

' A(p R = v 2 p R + 2uvipR - \(R)ipR, in B R (0), 

< Aip R = 2uv<pR + v 2 ipR-X{R)ip R , mB R (0), (3.2) 

k <PR = i>R = on dB R (0). 

By possibly replacing (pr,ipr) with — we can assume ifR > and ipR < 0. After a 

normalization, we also assume 

|^(0)| + |^(0)| = 1. (3.3) 
X(R) is decreasing in R, thus uniformly bounded as R — > +00. Let 

A:= Urn X(R). 

R^+oo 

The equation for <pr and — ipR (both of them are nonnegative functions) forms a cooperative 
system, thus by the Harnack inequality ([3] or [9]), pr and ipR are uniformly bounded on any 
compact set of M. 2 . By letting R — > +00, we can obtain a converging subsequence and the limit 
(ip, ip) satisfies (3.1). 

We also have ip > and ip < by passing to the limit. Hence 

-Aip+ {v 2 - X)<p > 0. 

Applying the strong maximum principle, either ip > strictly or p = 0. If p> = 0, substituting 
this into the first equation in (3.1), we obtain ip = 0. This contradicts the normalization condition 
(3.3). Thus, it holds true that tp > and similarly ip < 0. □ 

Fix a unit vector £. Differentiating the equation (1.10) yields the following equation for (u^,v^) 

{Auc = v 2 uc + 2uvve, 
Av^ = 2uvu£ + v v^. 

Let 

Wi = — ,W 2 = -j-. 

p Ip 

Direct calculations using (3.1) and (3.4) show 



{div(p 2 Vwi) = 2uvtpip(w2 — wi) + Xp 2 w±, 
div(p 2 ViV2) = 2uv(pip(wi — W2) + Xip 2 W2- 



11 



For any r\ £ Cq°(M?), testing these two equations with wirf and W217 2 respectively, we obtain 

— J ip 2 \Vwi\ 2 r] 2 — 2(p 2 wir]V w\\I r] = J 2uv(ptp(w2 — wi)wir] 2 + A (p 2 wir] 2 , 

— J ■lp 2 \Vw2\ 2 r] 2 — 2"4) 2 W2V1V W2^I t] = J 2uV(pif)(wi — W2)W2T] 2 + \ljj 2 W2t] 2 . 

Adding these two and applying the Cauchy-Schwarz inequality, we infer that 

J ^ 2 \Vwi\ V + t/> 2 |V™ 2 | V < 16 j <p 2 w 2 IVtjI 2 + i) 2 wl\Vi]\ 2 < 16 J (uf + v|)|V? ? | 2 . (3.5) 

Here we have taken away the positive term in the right hand side and used the fact that 
2uvtpijj(w2 — wi)wir] 2 + 2uvipip(wi — W2)w2T] 2 = —2uvipip(wi — W2) 2 rj 2 > 0, 

because 99 > and tp < 0. 

On the other hand, testing the equation Au > with urj 2 (77 as above) and integrating by 
parts, we get 



J |Vu|V < 16 J u 2 |Vr/p 



The same estimate also holds for v. For any r > 0, take 77 = 1 in B r (0), 77 = outside i?2r(0) and 
|Vr?| < 2/r. By the linear growth of u and v, we obtain a constant C such that 



\Vu\ 2 + \Vv\ 2 < Cr 2 . 



B r (0) 



(3.6) 



Now for any R > 0, in (3.5), we take 7/ to be 

7]{z) = 



h 
0, 



x € 5^(0), 



log 



^ .e%(0)\B B (0). 



With this 77, we infer from (3.5) 

/ tp 2 \Vwi\ 2 + i> 2 \Vw 2 \' 



< 



< 



'Bfl(O) 

c 
c 

(ioi^F 
c 



|Vd 2 + IVtjI 2 ) 



B b 2(0)\Bh(0) 

iJ 2 



I r~ 2 {! \Vu\ 2 + \Vv\ 2 )dr 

JR JdBr(0) 



f R r -^A. f |Vn| 2 + \Vv\ 2 )dr 
Jr dr J Br{0) 
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c 



,.Jr- 2 [ |Vn| 2 + \Vv\ 2 )\f + 2 f R r" 3 ( / |Vn| 2 + \Vv\ 2 )dr\ 

[VogK) JdB r (0) JR JB r (0) 

C 



< 



logR 

By letting R — > +00, we see Vw\ = and Vw2 = in M 2 . Thus, there is a constant c such that 

(ut,vt) = c(<p,if>). 

Because £ is an arbitrary unit vector, from this we actually know that after changing the coordi- 
nates suitably, 

u y = 0, v y = in M. 2 . 
That is, u and v depend on x only and they are one dimensional. 

4 Existence in bounded balls 

In this section we first construct a solution (u, v) to the problem 

( Au = uv 2 in flfl(0), 
1 Aw = to 2 in B R (0), 



(4.1) 



satisfying the boundary condition 

u = u = $ _ on &B fl (0) C K 2 . (4.2) 

More precisely, we prove 
Theorem 4.1. There exists a solution (ur,vr) to problem (4.1), satisfying 

1. Ufj — vr > in > 0} and — vr < in {<J> < 0}; 

2. Ur > $ + and v R > <E>~; 

3. Vi = 1, • • • , d, u R {Tiz) = v R (z); 
4- Vr € (0, 

, ^ B ,.(0)l V ^l 2 + l V ^l 2 + ^l , 
N(r;u R ,v R ) := ^ 2 2 < a. 

JdB r {0) U R + u il 



Proof. Let us denote C H 1 {Br{0)) 2 the set of pairs satisfying the boundary condition (4.2), 
together with conditions (1,2,3) of the statement of the Theorem (with the strict inequality < 
replaced by <, and so now U is a closed set). 
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The desired solution will be a minimizer of the energy functional 

E R (u,v):= / |Vu| 2 + |Vu| 2 + 

JBntO) 



2 2 
U V 



over U. Existence of at least one minimizer follows easily from the direct method of the Calculus 
of Variations. To prove that the minimizer also satisfies equation (4.1), we use the heat flow 
method. More precisely, we consider the following parabolic problem 



U t -AU = -UV 2 , in [0, +oo) x B R (0), 
V t -AV = -VU 2 , in [0, +oo) x B R (0), 



(4-3) 



with the boundary conditions U = & + and V = on (0, +oo) x dB R (0) and initial conditions 
inZ^. 

By the standard parabolic theory, there exists a unique local solution (U,V). Then by the 
maximum principle, < U < sup BR ^ <I> + , < V < sup Bfl ( ) <3?~, hence the solution can be 
extended to a global one, for all t £ (0, +oo). By noting the energy inequality 

and the fact that En > 0, standard parabolic theory implies that for any sequence ij — > +oo, 
there exists a subsequence of U such that (U(ti), V(ti)) converges to a solution (n, v) of (4.1). 

Next we show that IA is positively invariant by the parabolic flow. First of all, by the symmetry 
of initial and boundary data, (V(t, Tiz), U(t, Tiz)) is also a solution to the problem (4.3). By the 
uniqueness of solutions to the parabolic system (4.3), (U, V) inherits the symmetry of ( < l >+ ,$~). 
That is, for all t S [0, +co) and i = 1, • • • , d, 

U(t,z) = V(t,T iZ ). 

This implies 

U-V = on {$ = 0}. 
Thus, in the open set Dr := Br(0) n {3> > 0}, we have, for any initial datum (no, vq) € U, 

' (U - V) t - A(C7 -V) = UV{U - V), in [0, +oo) x D R {0), 
< U-V>0, on [0,+oo) x dD R (0), (4.5) 
k U-V > 0, on {0} x D R (0). 

The strong maximum principle implies U — V > in (0, +oo) x D R (0). By letting t — > +oo, we 
obtain that the limit satisfies 

u - v > in D R (0). (4.6) 
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(u, v) also has the symmetry, Vi = 1, • • • , d 

u(Tiz) = v(z). 



Similar to (4.5), noting (4.6), we have 

f -A(«-t;)>0, in D R (0), 
| u - v = $ + , on dD R (0). 

Comparing with <J> + on Dr(0), we obtain 

u - v > > 0, in D R (0). (4.8) 

Because u > and v > in .Br(O), we in fact have 

u > in J3fl(0). (4.9) 

In conclusion, (ti,f) satisfies conditions (1,2,3) in the statement of the theorem. 

Let (ur, vr) be a minimizer of Er over Now we consider the parabolic equation (4.3) with 
the initial condition 

U(x, t) = u R {x),V(x, t) = v R (x). (4.10) 

By (4.4), we deduce that 

Er(ur,vr) < E R {U,V) < Er(ur,vr) 

and hence (U(x,t),V(x,t)) = (ur(x),vr(x)) for all t > 0. By the arguments above, we see that 
(ur,vr) satisfies (4.1)and conditions (1,2,3) in the statement of the theorem. 

In order to prove (4), we firstly note that, as (ur, Vr) minimizes the energy and ( < & + , ( I >_ ) € U, 
there holds 

I \vu R \ 2 + \vv r \ 2 + 44< f |v$| 2 . 

Jb r (o) Jb r (o) 

Now by the Almgren monotonicity formula (Proposition 5.2 below) and the boundary conditions, 
Vr € (0, R), we derive 

RJ B m |V$| 2 
N(r; u R , v R ) < N(R; u R , v R ) < ( BrW ,^, 9 = d. 



JdB R (0) l^l 2 



>dB R (0) 

This completes the proof of Theorem 4.1. 



□ 



Let us now turn to the system with many components. In a similar way we shall prove the 
existence on bounded sets. Let d be an integer or a half-integer and 2d = hk be a multiple of the 
number of components k, and G denote the rotation of order 2d. Take the fundamental domain 
F of the rotations group of degree 2d, that is F = {z £ C : 9 = arg(z) € (— n/2d, ir/2d)}. 
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= r d cos(de) if z G U^G^F), (4n) 
1 otherwise in C. 

Note that ^(z) is positive whenever it is not zero. Next we construct a solution (u±, . . . ,u k ) to 
the system 

k 

A Ui = Ul u h in B R (0),i = l,...,k (4.12) 

satisfying the symmetry and boundary condition (here z is the complex conjugate of z) 

Ui{z) = Ui(G h z), on B R (0) ,i = l,...,k, 

Ui{z) = Ui+i(Gz), on B R (0),i = l,...,k, ^ ^ 

u k+2 -i{z) = Ui(z), on B R (0) ,i = l,...,k, 

u k+1 (z) = ui(z), on B R (0), 

Ui +1 (z) = V(G i (z)), on dB R (0) , i = 0, — , fe — 1. (4.14) 
More precisely, we prove the following. 

Theorem 4.2. For every R > 0, i/iere exists a solution (u± jR , . . . ,U}-,r) to the system (4.12) with 
symmetries (4.13) and boundary conditions (4.14), satisfying, 

r Ib <m \^ u i,R\ 2 + Ei<i u iR u2 jR 
N(r) := JBA0) — ' z < d, Vr G (0,12). 

JdB r (o) 2_a 

Proof. Let us denote by U C H 1 (B R (0)) k the set of pairs satisfying the symmetry and boundary 
condition (4.13), (4.14). The desired solution will be the minimizer of the energy functional 



/ J2 \^Ui,n\ 2 + V] uIruIr 



over U. Once more, to deal with the constraints, we may take advantage of the positive invariance 
of the associated heat flow: 

-^-AUi = -Ui^2uf, in [0,+ac)xB R (0), (415) 

which can be solved under conditions (1.20), (4.14) and initial conditions in U. Thus, the min- 
imizer of the energy {u\ tR , . . . , u k , R ) solves the differential system. In addition, using the test 
function (^i, . . . , Vf^), where V&j = ^ o G*" 1 , i = 1, . . . , k, we have 



k 

Jb r (o) V t^j j Br(o) 
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Now by the Almgren monotonicity formula below (Proposition 5.2) and the boundary conditions, 
we get 

R Jo m |Vlf 
N(r) < N(R) < JBr{0) = d, Vr G (0,1*). 

JdB R (0) 1*1 

□ 

In order to conclude the proof of Theorems 4.1 and 4.2, we need to find upper and lower 
bounds for the solutions, uniform with respect to R on bounded subsets of C. That is, we will 
prove that for any r > 0, there exists positive constants < c(r) < C(r) (independent of R) such 
that 

c(r) < sup ur < C(r). (4.16) 

Br(0) 

Once we have this estimate, then by letting R — > +oo, a subsequence of (ur, vr) will converge 
to a solution (u,v) of problem (1.10), uniformly on any compact set of M 2 . It is easily seen that 
properties (1), (2), (3) and (4) in Theorem 4.1 can be derived by passing to the limit, and we 
obtain the main results stated in Theorem 1.3 and 1.6. It then remains to establish the bound 
(4.16). In the next section, we shall obtain this estimate by using the monotonicity formula. 



5 Monotonicity formula 

Let us start by stating some monotonicity formulae for solutions to (1.10), for any dimension 
n > 2. The first two are well-known and we include them here for completeness. But we will also 
require some refinements. 

Proposition 5.1. For r > and x S W 1 , 

E{r)=r^l £|Vu 4 | 2 + ;>> 2 u 2 

JB r (x) 1 i<i 

is nondecreasing in r. 

For a proof, see [12]. The next statement is an Almgren-type monotonicity formula with 
remainder. 

Proposition 5.2. For r > and x S W 1 , let us define 

JdBr(x) I 

Then 

iV(nx):=§l 

H[r) 
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is nondecreasing in r. In addition there holds 



I 

Jo 



2 Sb s ^i<j 



k u 2 u 2 



<3 * 3 



IdB, £l 



k „,2 



ds < N(r) . 



(5.1) 



Proof. For simplicity, take x to be the origin and let k = 2. We have 



H(r)=r 1 - n [ u 2 +v 2 , E{r) = r 2 ~ n I \Vu\ 2 + | Vv\ 2 + ttV 
J dB r J B r 

Then, direct calculations show that 

-f-tf(r) = 2r l ~ n f \Vu\ 2 + \Vv\ 2 + 2u 2 v 2 . 
dr J Br 

By the proof of Proposition 5.1, we have 

-E(r) = 2r 2 ~ n f [u 2 . + v 2 .} + 2r x ~ n j 

J dB r B r 



d_ 

dr 



2 2 

u v . 



(5.2) 



(5.3) 



With these two identities, we obtain 

±E = H[2r*~ n f dBr (u 2 + v 2 ) + 2r^» u 2 v 2 } - E[2r^ nn r + iny 
drF (rj H 2 



> 



^ 2n I 9Br {u 2 + v 2 ) j dBr (u 2 + v 2 )- 2r 3 " 2 " J dBr uu T + 



vr r 



H 2 



-+ 



1-n C U 2 V 2 



+ 



2r l-nj i 



H 



> 



2r^J i 



u 2 v 2 



H 



Here we have used the following inequality 

E(r) < / |Vu| 2 + |Vt;| 2 + 2u 2 v 2 = / uu r + 

J B r J 0B r 



vr r 



Hence this yields monotonicity of the Almgren quotient. In addition, by integrating the above 
inequality we obtain 



2„,2 



IdB s u2 + V ' 



: ds < N(r) . 



□ 



If x = 0, we simply denote N(r; x) as N(r). Assuming an upper bound on N(r), we establish 
a doubling property by the Almgren monotonicity formula. 
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Proposition 5.3. Let R > 1 and let (u±, . . . , Uk) be a solution of (1.19) on Br. If N(R) < d, 
then for any 1 < r\ < r<i < R 



H ( r 2) < d r 2 

H{n) ~ 



2,1 



2d' 



(5.4) 



Proof. For simplicity of notation, we expose the proof for the case of two components. By direct 
calculation using (5.2), we obtain 



d i 

dr 



u 2 + v 2 ) 



dB r (0) 



2 f Br \Vu\ 2 + \Vv\ 2 + 2 



u 2 v 2 



JdB r (0) U<2 + y2 



< 



2N(r) , 2/j 



2„,2 



+ 



B, 



U V 



2d 
< — + 



JdB r (0) 



n 2 + v 2 



JdB r (0) 



u 2 + v 2 



Thanks to (5.1), by integrating, we find that, if T\ < r<i < 2rg then 



H{r 2 ) 



r 2d 
— Jld 



(5.5) 



□ 



An immediate consequence of Proposition 5.3 is the lower bound on bounded sets for the 
solutions found in Theorems 4.1 and 4.2. 

Proposition 5.4. Ler (ui t R, . . . , u^ r) be a family of solutions to (1.19) such that N(R) < d and 
H(R) = CR 2d . Then, for every fixed r < R, there holds 

H(r) > Ce- d r 2d . 



Another byproduct of the monotonicity formula with the remainder (5.1) is the existence of 
the limit of H(r)/r 2d . 

Corollary 5.5. Let R > 1 and let (m, . . . , Uk) be a solution of (1.19) on C such that limr-^+oo N{r) < 
d, then there exists 

lim —kJ- < +oo . (5.6) 



r— >+oo r 



2d 



Now we prove the optimal lower bound on the growth of the solution. To this aim, we need a 
fine estimate on the asymptotics of the lowest eigenvalue as the competition term diverges. The 
following result is an extension of Theorem 1.6 in [8], where the estimate was proved in case of 
two components. 
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Theorem 5.6. Let d be a fixed integer and let us consider 

c2tt d d 



C(d, A) = min < 



J i i<3 



2 2 



Jo^ Si ^ 2 = 1) u i+1 (x) = Ui(x - 2n/d), 

ui(-x) = Ui(x) ,U d+1 = U\ 



(5.7) 



Then, there exists a constant C such that for all A > 1 we have 

d 2 -CA~ 1/4 < C(d,A) < d 2 . 
Proof. Any minimizer (tti a? • • • i u d,h) solves the system of ordinary differential equations 



II; 



Uj — \Ui 



L j . . . , W,, 



d, 



(5. 



(5.9) 



together with the associated energy conservation law 



^2(u' t ) 2 + An 2 - A ^ u 2 u 2 = h . 



(5.10) 



Kj 



Note that the Lagrange multiplier satisfies 

r-2-K d d 



A 



r x: ki 2 + 2A £ ^ = £( d> a) + r a ^ 

jo r~ JO 



As A — >• oo, we see convergence of the eigenvalues A ~ £((f, A) — > <i 2 , together with the energies 
h — > 2d 2 . Moreover, the solutions remain bounded in Lipschitz norm and converge in Sobolev 
and Holder spaces (see [8] for more details). Now, let us focus on the interval / = (a, a + 2n/d) 
where the i-th component is active. The symmetry constraints imply 

Ui-i(a) = Ui(a) ,ii;_i(a) = -u'^a) , 

Ui + \(a + 2-k jd) = Ui(a + 2n/d) , u i+1 (a + 2-k /d) = —u^a + 2ir/d) 

We observe that there is interaction only with the two prime neighboring components, while the 
others are exponentially small (in A) on I. Close to the endpoint a, the component m is increasing 
and convex, while U{-\ is decreasing and again convex. Similarly to [8] we have that 



u 



(a) = Ui-i(a) ~ KA" 1 / 4 X(a) 



H = (h + K)/2 



(5.11) 



Hence, in a right neighborhood of a, there holds Ui(x) > Uj(a), and therefore, as u i _ 1 > 
Au 2 (a)uj_i, from the initial value problem (5.11) we infer 



U:. 



-l(x) < C Ui (a)e- Al/2u ^ a){x - a) ,Vx € [a, b]. 
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On the other hand, on the same interval we have 

Ui(x) < Ui(a) + C (x — a) ,Vx G [a, b] . 
(here and below C denotes a constant independent of A). Consequently, there holds 

A J^uf^u] + uf_ lUi + m-mf < CA'^Uiia)- 1 ~ CA" 1 / 4 . (5.12) 

In particular, this yields 

C(d,A) > X-CA- 1 / 4 . 



In order to estimate A, let us consider m 



±1 3 



(5.13) 

Then, as Ui{a) = Uj_i(a) and 

Ui{a + 2-7T / 'd) = Ui + \(a + 2'K / d), Ui G H l {I). By testing the differential equation for Ui — Ej=i±i u j 
with Ui on / we find 

jf|S;[ a <Ajf[Tfc| a + C7A- 1 /*, 

where in the last term we have majorized all the integrals of mixed fourth order monomials with 
(5.12). As |/| = 2n/d, using Poincare inequality and (5.13) we obtain the desired estimate on 
C(d,A). □ 

We are now ready to apply the estimate from below on £ to derive a lower bound on the 
energy growth. We recall that there holds 



E(r):=[ £|Vu,| 2 + 2]Tu^= / J> 

JB r {x) ~ JdB r (x) 



duj 
dr 



Proposition 5.7. Let (u± r, . . . , r) be a solution of (1.19) having the symmetries (1.20) on 
Br . There exists a constant C (independent of R) such that for all 1 < t\ < r<i < R there holds 



E(r 2 ) 

E( n ) 



> C- 



,:2d 



.2d 



(5.14) 



Proof. Let us compute, 



d . ( A 2d , f dBr(x) Ei |V^| 2 + 2E i <i^i 



— log ( r~ 2d E(r 
dr 



— + 

r 



IdB r (x) El 



Ui- 



2d 
+ 

r 



fdB r {x) El 



duj 
dr 



+ 



duj 
dr 



duj 



+ 2r 2 E l<J uW j 



IdB r (x) El 



U i 



duj 
dr 
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2d 



Jo z2i Ui dr 



Now we use Theorem 5.6 and we continue the chain of inequalities: 

2 r(A 



t-2-K ( 9ui\ £(d,2r ) „ 27r „ fe 2 
2d Jo ^{dr J + ~ r2 J o ^1 Ui 



-log r^E(r) > + - 

dr \ w 7 r r 27r ^fc 

Jo Z^i M * g r 

a _ M -vro a __^ ( , 15) 

where in the last line we have used Holder inequality. By integration we easily obtain the assertion. 

□ 

A direct consequence of the above inequalities is the non vanishing of the quotient E/r 2d : 

Corollary 5.8. Let R > 1 and let (u\, . . . , it*.) be a solution of (1.19) on C satisfying 1.20: then 
there exists ^ 

lim = b£ (0,+oo] . (5.16) 

If, in addition, lim r _>._|_ 00 N(r) < d, then we have that b < +oo and 

E(r) 

lim N(r) = d, and lim = b . (5.17) 

r— H-oo i — ^+oo T"'™ 

Proof. Note that (5.16) is a straightforward consequence of the monotonicity formula (5.15). To 
prove (5.17), we first notice that 

lim ^1 = lim N(r)*£. 



So the limit of E(r)/r 2d exists finite. Now we use (5.1) 



2„,2 



JdB s Si U i 

and we infer 



ds < +oo 



,2„,2 



liminf r/ ^ E ^^ =0. 



fc 2 



22 



Next, using Corollary 5.5 we can compute 



lim inf UTLdA _ lim M J^^ _ o, 

r— >-+oo r ia r— >+oo H(r) r za 



and finally 



iiminfE(r) E(r) _ 

r— ¥+00 f 



Was the limit of N(r) strictly less that d, the growth of H(r) would be in contradiction with that 
of E(r). □ 

Now we can combine the upper and lower estimates to obtain convergence of the approximating 
solutions on compact sets and complete the proof of Theorems 1.6 

Proof of Theorem 1.6. Let (ui r, . . . , Uk r) be a family of solutions to (1.19) such that Nr(R) < d 
and Hr(R) = CR 2d . Since Hr(R) = CR 2d , then, by Proposition 5.3 we deduce that, for every 
fixed 1 < r < R, there holds 

H R (r) > Ce~ d r 2d . 
Assume first that there holds a uniform bound for some r > 1, 

H R {r) < C . (5.18) 

Then HR{r) and Er{t) are uniformly bounded on R. This implies a uniform bound on the 
i? 1 (-B r ) norm. As the components are subharmonic, standard elliptic estimates (Harnack inequal- 
ity) yield actually a C 2 bound on B r / 2 , which is independent on R. Note that, by Proposition 5.4, 
HR(r) is bounded away from zero, so the weak limit cannot be zero. By the doubling Property 
5.3 the uniform bound on Hr(v2) < Cr 2d holds for every r 2 € M larger than r. Thus, a diag- 
onal procedure yields existence of a nontrivial limit solution of the differential system, defined 
on the whole of C. It is worthwhile noticing that this solution inherits all the symmetries of the 
approximating solutions together with the upper bound on the Almgren's quotient. Finally, from 
Corollary 5.5 and 5.8 infer the limit 

H(r) 1 E(r) b 

lim — £p = lim "TTTT' lim -4/ = -j e (0, +00) . (5.19) 
r r—t+oo Jy If j r— >+oo r d 



-00 



Let us now show that Hr(t) is uniformly bounded with respect to R for fixed r. We argue by 
contradiction and assume that, for a sequence R n — > +00, there holds 

lim H Rn {r) = +00 . (5.20) 

n— >+oo 

Denote Ui tTl = U{ t R n and H n , E n , N n the corresponding functions. Note that, as E n is bounded, 
we must have N n (r) — > 0. For each n, let A n G (0,r) such that 
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(such A n exist right because of (5.20)) and scale 

n z) , \z\ < R n /X n . 

Note that the (ui, n )i still solve system (1.19) on the disk B(0, R n /X n ) and enjoy all the symmetries 
(1.20). Let us denote H n , E n , N n the corresponding quantities. We have 

H n (l) = XnH n (X n ) = 1, 
E n (l) = X 2 n E n (X n ) -> 
N n (l) = N n (X n ) -> 

In addition there holds N n (s) < d for s < R n /X n . By the compactness argument exposed above, 
we can extract a subsequence converging in the compact-open topology of C 2 to a nontrivial 
symmetric solution of (1.19) with Almgren quotient vanishing constantly. Thus, such solution 
should be a nonzero constant in each component, but constant solution are not compatible with 
the system of PDE's (1.19) . □ 



6 Asymptotics at infinity 

We now come to the proof of Theorem 1.4. Note that by Proposition 5.3, the condition on N(r) 
implies that u and v have a polynomial growth. (In fact, with more effort we can show the reverse 
also holds. Namely, if u and v have polynomial growth, then N(r) approaches a positive integer 
as r — > +oo. We leave out the proof.) 

Recall the blow down sequence is defined by 

(u R (x),v R (x)) := (j^ u ( Rx )>J^ v ( Rx ))' 



where L(R) is chosen so that 



4 + 4= / <*> 2 - (6.i) 

aBi(O) JdBtiO) 



For the solutions in Theorem 1.3, by (5.19), we have 

L{R) ~ R d . (6.2) 
We will now analyze the limit of (ur,vr) as R — > +oo. 

Because for any r E (0, +oo), N(r) < d, (u,v) satisfies Proposition 5.3 for any r £ (l,+oo). 

After rescaling, we see that Proposition 5.3 holds for (ur, vr) as well. Hence, there exists a 
constant C > 0, such that for any R and r € (1, +oo), 

4 + 4<CeV. (6.3) 

8B r (0) 
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Next, (ur, vr) satisfies the equation 

Au R = L(R) 2 R 2 urv 2 r , 



Av R = L(R) R vru r , (6.4) 



Ur, v r > in R 2 . 



Here we need to observe that, by (6.2) 



lim L(R) 2 R 2 = +oo. 

By (6.3), as R — > +oo, ur and vr are uniformly bounded on any compact set of ]R 2 . Then by 
the main result in [16], [24] and [26], there is a harmonic function \E' defined in IR 2 , such that (a 
subsequence of) {ur, vr) — > (\& + , in H 1 and in Holder spaces on any compact set of M 2 . By 
(6.1), 



/ * 2= / 

JdBi (0) Jdl 



$2 



'afii(O) J9Bi(0) 
so ^ is nonzero. Because L(R) — > +oo, ur(0) and vr(0) goes to 0, hence 

9(0) = 0. (6.5) 

After rescaling in Proposition 5.2, we obtain a corresponding monotonicity formula for (ur, vr), 

rf Br{0) \VuR\ 2 + \Vv R \ 2 + L(R) 2 R 2 u 2 R v 2 R 
N(r;uR,VR) := — 2 2 = N(Rr) 

JdB r (0) U R + V R 

is nondecreasing in r. By (4) in Theorem 1.3 and from Corollary 5.8, 

N(r;u R ,VR) <d= lim N(r;u R ,v R ) ,VrG(0,+oo). (6.6) 

r— >+oo 

In [16], it's also proved that (ur,vr) — > (^ r+ ,^ r_ ) in Hj oc and for any r < +oo, 

lim / L(R) 2 R 2 u 2 R v R = 0. 

After letting i? — > +oo in (6.6), we get 

r f [V^l 2 

N{r;V):= f r(0) — = lim N(r;u R ,v R ) = lim N{Rr) = d. (6.7) 

W(0) ^ ^ + °° ^ + °° 

In particular, iV(r; \&) is a constant for all r £ (0,+oo). So is a homogeneous polynomial 
of degree d. Actually the number <i is the vanishing order of at 0, which must therefore be 
a positive integer. Now it remains to prove that \& = $: this is easily done by exploiting the 
symmetry conditions on 9 (point (3) of Theorem 1.3). 
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